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Abstract 

We consider the orthogonal polynomial {P n (z)} with respect to the measure 

2 — a| 2JVc e _JV ' z ' dA(z) 

over the whole complex plane. We obtain the strong asymptotic of the orthogonal polynomials 
in the complex plane and the location of their zeroes in a scaling limit where n grows to infinity 
with N. 

The asymptotics are described in terms of three (probability) measures associated with the 
problem. The first measure is the limit of the counting measure of zeroes of the polynomials, 
which is captured by the g-function much in the spirit of ordinary orthogonal polynomials 
on the real line. The second measure is the equilibrium measure that minimizes a certain 
logarithmic potential energy, supported on a region K of the complex plane. The third measure 
is the harmonic measure of K c with a pole at oo. This appears as the limit of the probability 
measure given (up to the normalization constant) by the squared modulus of the n-th orthogonal 
polynomial times the orthogonality measure, i.e. |P n (z)| 2 |,z — a| 2JVc e _JV ' z ' dA(z). 

The compact region K that is the support of the second measure undergoes a topological 
transition under the variation of the parameter t = n/N; in a double scaling limit near the 
critical point given by t c = a(a + we observe the Hastings-McLeod solution to Painleve 

II in the asymptotics of the orthogonal polynomials. 
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1 Introduction 

In this paper, we study orthogonal polynomials with respect to a weight function that is supported 
on the whole complex plane. Let P n {z) = P n .N( z ) be the monic orthogonal polynomial such that 

/ P n , N (z) P miN (z) e- NQ ^dA(z) = h n , N 5 nm (rc,m = 0,l,...), (1.1) 
Jc 

where N is a positive real parameter, Q : C — > R is called the external potential, dA(z) is the 
standard area measure on C, and h n = h n ^ is the (positive) norming constant. 

Our principal motivation is the connection between these orthogonal polynomials and so-called 
"random normal matrices". (See [24] for a complete description of this connection.) Briefly, one 
considers the space of n x n normal matrices (those that commute with their adjoint) equipped 
with a probability measure of the form 

1 e -NT*[Q(M,M')] dM , (1.2) 

where dM represents the measure induced on normal matrices by Lebesgue measure on the space 
C n of all n x n complex matrices. One of the fundamental quests within random matrix theory 
is to understand the statistical behavior of the eigenvalues in the limit when n (and N) tend to 
infinity. The amazing fact is that most natural statistical quantities concerning eigenvalues can be 
expressed in terms of the associated orthogonal polynomials. In particular, the averaged density of 
eigenvalues is precisely 

a n M^-e- NQ ^J2 lP f^. (1.3) 

3=0 J ' 

Using this fact, once one has complete asymptotic control of the orthogonal polynomials, it is an 
interesting challenge to compute many different statistical properties of the eigenvalues in a concrete 
fashion as n and N tend to oo. 

There is only a handful of external potentials Q for which these polynomials can be explicitly 
computed, and the average density can be analyzed asymptotically, as n and iV tend to oo. Taking 
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the simplest case, Q(z) = \z\ 2 , one obtains the hj^'s and the averaged density of eigenvalues 0"„,jv 
in terms of the Gamma function. For this choice of Q, the associated matrix model is called the 
Ginibre ensemble, and has been studied to great effect via the simple orthogonal polynomials z n 
(see for instance [24, 17] and the appendix of [27]). In fact, for any radially symmetric potential 
Q(z) = Q(\z\), the orthogonal polynomials are the monomials z n . 

A harmonic deformation of the Gaussian case, i.e. Q(z) = \z\ 2 + (harmonic), is particularly 
interesting (because of its relation to a growth model called the Hele-Shaw flow [26]). The simplest 
such deformation of the Gaussian normal matrix model uses the potential Q = \z\ 2 + t (z 2 + z 2 ) . 
The associated orthogonal polynomials are now Hermite polynomials [15], for which the asymptotics 
are well known. The case with cubic deformation was recently studied in [8]. 

Concerning the limiting density of eigenvalues, there are other approaches [16, 2], which connect 
to a classical energy problem from potential theory. 1 

Let us scale the real parameter TV with n so that lim„_ s . 00 n/N — t for some fixed positive number 
t. The general fact [16, 2] is that, under some basic assumptions on the potential Q, the limiting 
density of eigenvalues converges (in the sense of measures) to the unique probability measure [i* in 
the plane which minimizes the functional C(fj,): 

£(A0 - * / Q{z) M*) + t 2 (T log r-^—rd^z) dfi(w). (1.4) 

This functional is the well-known Coulomb energy functional in M 2 , and has a well-developed 
classical variational theory [25] . The measure fj,* is equivalently characterized in terms of the Eulcr- 
Lagrange variational conditions, which can be summarized as follows: fi* is said to satisfy the 
Euler-Lagrange conditions if there is a constant £ 2 d so that 

Q{z)-2t [ \og\z-w\ dfi* (w) + £ 2 d > 0, (1.5) 
Jc 

for all z in C (up to a set of capacity zero), with equality on the support of the measure fi* (again up 
to a set of capacity zero). The left hand side of the above is denoted by U2D, the effective potential 
of the Coulomb system, and of which the smooth extension, Li, (defined at (2.24)) will show up in 
Corollary 1.7 in a conspicuous way. 

In this paper we study one of the simplest cases, given by 

Q(z) := \z\ 2 + 2clog - c>0, a>0. (1.6) 
\z — a\ 

The case of a 6 C \ {0} can be reduced to the above by a rotation of the coordinate. We will assume 
a^O (for a = 0, one has P n ,N{z) — z n ). 

This case has been studied in [3] where the corresponding equilibrium measure is described. In 
Section 2 we define the set K C C and prove that it is the support of the equilibrium measure 
(Proposition 2.4 and Proposition 2.6). The region K is a function of the parameter t such that the 
area of K is given by irt. As Figure 1 indicates, K undergoes a topological transition at the critical 
time t c given by 

t c := a(a + 2-y/c). (1.7) 
With respect to this, we separate our analysis into the three cases: 



In random Hermitian matrices, these are analogous to the saddle point approach from the physics community 
(see, for example [9]), which is put on a solid mathematical footing originally by Johansson [20]. 
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• the pre-critical case, < t < t c : The support set K is a simply connected region with 
analytic boundary. The left-hand figure of Figure 1 shows a representative picture for this 
case, and Lemma 2.5 describes K explicitly using a conformal map. 

• the post-critical case, t > t c : K consists of a region between two circles, a larger circle 
centered at the origin and a smaller one centered at a. The right-hand figure of Figure 1 
shows a representative picture, and Proposition 2.4 describes the domain K. 

• the critical case, t ~ t c : where we consider a scaling limit to the criticality such that 

t-t c = (D(N- 2/3 ). (1.8) 



Figure 1: K (shaded regions), B (lines) and (3 and/or j3 (dots). The pre-critical (left), critical 
(middle), post-critical (right) cases. 



Summary of results: We consider the following version of the scaling limit, given by 

N = (n — r)/t, n oo, for fixed constants, t £ ]R + and r e Z. (1-9) 

We express our asymptotic results in terms of geometric quantities related to K (we believe that, 
in such form, the asymptotic expression holds for more general Q). Let us define ~P{K) to be the 
polynomial convex hull of K (which in the present cases may be defined as the set whose complement 
is precisely the connected component of K c containing oo). It will turn out that K is either simply 
connected (t < t c ), or multiply connected (t > t c ) as shown in Figure 1, and so 

P(K) = { K \ , 1 P re " CTitical case - (L10 ) 

I {z : \z\ < \/t + cj, critical and post-critical cases. 

Below, we describe some quantities that arc determined by the geometry of K. 

- Let F be the conformal map F : C \ P(K) — >• C \ D (where D is the unit disk centered at 
the origin and C = C U {oo}) such that F(z) — z/p + 0(1) for some p > as z — > oo. The 
constant p is called the logarithmic capacity ofP(K). 
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- For the pre-critical case, the branch cut B (Definition 2.15) is a curve within K between two 
endpoints that we denote by (3 and B. For the post-critical and the critical cases, B (Definition 
2.19) is a Jordan curve passing through (3 £ R. We can therefore define Int(S) and Ext(S) 
such that oo € Ext(B). 

- We extend the domain of F: In the pre-critical case, we will show that F extends analytically 
to C \ B. In the post-critical well as the critical case, we show that F extends 
analytically to Ext(S). 

- Let g(z) be defined in Definition 3.5 (in terms of y that is defined for all three cases in Section 
2). For z $l K, it is exactly the complex logarithmic potential of the measure uniformly 
supported on K (see Lemma 3.7). For the critical and the post-critical cases, g{z) is given 
explicitly by (3.27). The complex valued function <j), which is purely imaginary on B, is defined 
in Definition 3.4. 

To state our results, we also note that Q,± is defined at Definition 3.8. (Also see Figure 7 and Figure 
8.) We further define a few constants: 



^■-f-^r/y b c :=a + V-c, (1.11) 

and expicit formulae for B and b (for the post-critical case) are given in (2.31). 

Theorem 1.1 (Pre-critical). For any fixed t < t c , let K be the region defined in Lemma 2.5. We 
have 

P n A^)^^ + 0{N- 1 ))^pW{z){pF{z)Ye m ^\ (1.12) 
uniformly over z in compact subsets of C\B. For z near B but away from B and B, we have 

P„,„ W = e'«-> (orfW ± e«<'» (iL)' + 0(W -,) , (LIS, 

uniformly over z in compact subsets of Q±\{[3, [3} (respectively for ±). The real parameters, k and 
a, are uniquely determined from the conformal map F by 

F(z) = - + — + -+0 (-5- , z^oo. 1.14 
p ap z \z 1 ) 

Lastly, to describe P n ,N near [3, let us define z(Q := [3 + (C^-iV)~ 2 / 3 C with 



C r -= ■ (1-15) 



a( f3-b)Vf3- [3 
W -o) 

where f3 and b are defined at (2. 12) (2. BO). We have 

Pn,N(z(0) - ( 4 . ( ^_ ) V4 (v.pYVnC 1 ^ (Ai(C) + O (iV-V3) ) e IC 3/2 e^^(0), (i.l 6 ) 

uniformly over ^ in compact subsets of C. The correct branch of the 4th-root is determined by 
knowing that the factor comes from the leading behavior of \J pF'(z) when z — > 8. The case for 
z — > B is obtained by the replacement: j3. 
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In (1.16) we write the asymptotic formula such that one can easily see that it matches the 
formula (1.12) as \((z)\ — > oo. 

Theorem 1.2 (Post-critical). Fort > t C! we have 



/ \ JVC 



uniformly over z in compact sets ofC\B. The constant t is defined in (3.28). For z near B but 
away from f3, we have 

p n {z) = { \ z - a J v v^Nz-( Z -f3) v w*-- (118) 



uniformly over z in compact subsets of Q± \ {/3}. 

To describe P n ^ near j3, we define z(£) := (3 + 7 1 iV~ 1 / 2 £ with 'ji defined at (1.11). We have 

PnA*(0) = e Naz+tm p r (j-(C) + O (V 1 / 2 )) , (1.19) 
uniformly over £ in compact subsets of C, where J- is the entire function given by 

1 f°° e s2/2 , fo, 3RC < 0. , nn . 

(This function is plotted in Figure 10.) 

Theorem 1.3 (Critical). In the critical case, i.e. t — t c = (D(N~ 2 / 3 ), we have 

\ / \ Nc 

+ o (*-"■))=• (^ 



uniformly over z in compact subsets of C \ B. Here s = c 1 / 6 a _1 / 3 (6 c )^ 2 / 3 A^ 2 / 3 (i — t c ) and u(s) := 
q'(s) 2 — s q(s) 2 — q(s) 4 with q(s) being the Hastings-McLeod solution of Painleve II equation. 
For z near B but away from ft, we have 



Pn(z) = { 



Nc 

' ~ 2N 1 / 3 (z-b c ) ' 2N 1 / 3 z r (z-b c ) 



±_\ c ( 1 _ _jM 3 ) + M^H^gj^ e^w + o (at- 2 / 3 ) ) , 



2NV 3 (z-b c )J 2N 1 / 3 {z-b c ) 

— — ' (L22) 

uniformly over z in compact subsets, respectively, of Q— \ {/3} and of \ {/?}. 
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Remark. We note that the leading asymptotic behavior of P n in (1.12) is conveniently expressed 
in terms of the geometric quantities, p, F and g. The asymptotic formula in Ext (6) for the post- 
critical and critical cases shown in (1.17) and (1.21) can also be expressed in terms of the geometric 
quantities if one uses pF(z) = z and (3.27). Then the leading terms take precisely the same 
form as the leading term in (1.12), i.e. z n (z/(z — a)) Nc = (pF(z)) r y/pF'(z) e Ntg ( z \ It has been 
conjectured [1] that the asymptotic behavior in K c should always take the latter form. 

Remark. One can confirm that the strong asymptotics in the critical case match those in the 
post-critical case in the limit s — > oo, using 

9(s)e<fC| ~2^?74 *"»->°°> and 7*~^vl' (1-23) 

where the former is obtained from the asymptotic behavior of the Hastings-McLeod solution (6.15), 
and the latter is by writing 7 (1.11) using (2.32). An argument for the appearance of the Painleve 
II equation in the critical case was explained in [22]. 

Remark. In the above theorems, though not immediately noticeable, the leading asymptotics 
(1.13)(1.18)(1.22) near B are analytic on B \ {(3 U 0) and they are exactly the sum of the two 
contributions from both sides of B by analytic continuations of the asymptotic formula on each side. 
We also note that the two leading terms for z E f2± are simply related by [tg(z)}± — [tg(z) + <j>(z)]^. 

Equipped with the above asymptotic results concerning the orthogonal polynomials, it is straight- 
forward but tedious to prove the following proposition. In this paper we do not give a complete 
proof of this proposition. Rather, we provide (at the end of Section 5) a sketch of the important 
steps involved in the proof. (See [21] for similar argument regarding the location of zeroes.) 

Proposition 1.4. The support of the counting measure of the zeroes of the polynomials P ny N 
outside of arbitrary small disk centered at the turning points z = f3,fi (pre- critical) and z = f3 
( critical/ po stcritical) , tends uniformly to the branchcut B for all three cases. Moreover, in the post- 
critical case, zeros are within 0{N~ 1 ) away from the curve defined implicitly by (5.24) that tends 
to B at the rate of OQxi N/N) (see Figures 11 and 12). 

It is possible to deduce from the proof of Theorem 1.3 an asymptotic description for P rit N which 
is valid near /? for the critical case. The formulae, however, are cumbersome in that they involve 
the solution of the Painleve II Riemann-Hilbert problem itself. Rather than highlighting this with 
a separate theorem, we provide instead the following proposition, establishing a uniform bound on 
P n ,N m t ne entire plane. 

Proposition 1.5. We have the following uniform bounds: 

0(iV _1 / 2 ), over compact subsets oflnt(B) for post-critical; 

n jvtg(^) 1 OjW -1 / 3 ), over compact subsets of Int(£>) for critical; (124) 

©(iV 1 / 6 ), on a neighborhood of '{/3 , /?} for pre- critical; 

v O(l), otherwise. 

The following theorem, that we prove in Section 7, gives the asymptotics of h n . 
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Theorem 1.6. For all t, we have 

h n ={l + 0(N-^))2^^ N p^e N ^, (1.25) 
where £20 is the real number given by £20 — 2t$lg(z) — Q{z) for z £ dP(K). 




Figure 2: The graph of 2 arctan(W(z)) for the pre-critical case; the right side is a view from the top. 
The parameters are a = 3.7619, c = 6.9168, t = 4.0557 (p = 2.1; a = 0.4; k = 1/2). Indicated are 
the boundary of K (thin green line) and the branch cut B. 



Harmonic measure from P n .N- The harmonic measure ui of a domain C C with a pole at 
00 (assume 00 e Q) is the probability measure supported on dft that satisfies, for any continuous 
function / on dft, the identity / (00) = f(z) dw(z) where / : CI — > C is the unique bounded 
harmonic extension of /. If d£l consists of piecewise smooth Jordan curves, then we have [18] 

Aw{z) = ^\F\z)\dI{z) (1.26) 

where F : fl — > D is a conformal map that preserves 00, and d£ is the length measure on dfl. 
We define the following probability density function: 

Pn (z):=±\P n (z)fe- N ^. (1.27) 

K 

This is but a single term in the average density (1.3) of the eigenvalues, representing "the density 
of an additional eigenvalue on top of n eigenvalues". In the scaling limit, an additional eigenvalue 
gives an infinitesimal increase to t by an amount of 1/N (recall that t ~ n/N). By considering the 
pertubation of t in the functional C{p) (or, equivalently, by taking the i-derivative of (1-5)), one can 
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f 



Figure 3: Plot of 2tt^tt/(2N) p n (z) for c = 1/6, a = 1, N = 30, n = 25, t = 0.83, i c R 
1.81649, Nt c — 54 ~ 55. The red dot is the location of z = a; The red line shows the height \F'(z) 
for z £ dK . 



deduce that the incremental growth of the measure, d(£/x*)/dt, is in terms of the harmonic measure 
of C \ ~P{K) with a pole at oo. This suggests that the measure p n dA converges to such harmonic 
measure. This fact was conjectured in physics literature [1] and proved in Theorem 7.7.3 of [2]. 

Below we obtain global uniform asymptotics for p n , which, not only shows that p n dA converges 
to the harmonic measure, but also shows in a very precise way how p n converges pointwise. We 
mention that the leading term in (1.27) was reported in the physics literature [1]. 

Corollary 1.7. We have p n (z) — 0{l)eT Nu ^ z ' uniformly in z £ C and, furthermore, 

Pn (z) = (i + o (at-i/3)) m^^E w > [z)]e -NU( Z) (L28) 

uniformly over z in compact subsets ofC\B for the pre-critical case and o/Ext(i3) for the critical 
and post-critical cases. Over compact subsets in Int(B), we have p n {z) — 0{N~ V ) q~ nu ( z ) for the 
post-critical case, and p n {z) — 0(N~ 2 / 3 ) qt nu ( z ) for the critical case. For the pre-critical case, 
over a neighborhood ofB, we have p n (z) = 0(N 1 ^ 3 ) e~ NU ( z ) . 

Proof. The proof follows from Proposition 1.5, the asymptotic formula (1.12) and Theorem 1.6. Let 
us only check it for z in compact subset of C \ 23 and of Ext (23) (the other regions can be similarly 
considered). Using the remark below Theorem 1.3 the leading term (1.12) applies to all values of t. 
Plugging in the asymptotics of P n (1.12) and h n (1-25) into p n (1-27) we obtain 

p n (z) = (l + (N-V>)) ^g^e-^) (1.29) 

The error bound that comes from the worst case: (1.21) in the critical case. The result follows from 
the identity U(z) — Q(z) — 2tdig(z) + £21) which can be seen from (2.2) and Lemma 3.7. □ 
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Corollary 1.7 says that, for the pre- and the post-critical cases, the probability density function 
p n is exponentially suppressed away from dP(K) where U is minimized (Lemma 2.18, Lemma 2.20). 
For the critical case, U is minimized on dK by Lemma 2.21 so p n may be peaked on dK n Int(£>) 
also. However, p n = 0{N~ 2 / 3 ')e~ NU ^ in Int(S), and this means only 0(N~ 2 ^ 3 ) amount of mass 
from the probability measure p n dA can be within Int(£>). For the pre-critical case, no mass can 
be accumulated at /3 or /3 in the scaling limit because 0{N 1 / 3 ') q~ nu ( z ) is exponentially suppressed 
there. To see how we obtain the approximation of the harmonic measure, let p be a point on dP(K) 
and n be a unimodular complex number whose direction in the complex plane is perpendicular to 
dP(K). Then we see from (2.6), the fact that U is the smooth extension of W2D to the interior, 
together with (1.28) that 

Pn (p + X n) ~ \ F \ +X ^ \F( P + X n)rJ™e-™* 2 (1.30) 

where we see the normal distribution approximation of the Dirac delta function of the boundary 
[x = 0) in the trailing term. This means that the density p n (p+xn) is localized around dP(K) with 
the typical width ~ y/l/N. Since ^\F'(p)\d£(p) (with d£ being the length measure on dP(K)) 
gives the harmonic measure and |-F(p)| = 1, we conclude that p n dA converges to the harmonic 
measure and moreover p n converges pointwise to the density of the harmonic measure (with respect 
to the area measure). 

Strategy of the proofs: In Section 2, we begin by defining the region K in terms of explicit 
formulae. We introduce the Schwarz function S(z) in terms of which we define y{z) that facilitates 
the construction of the ^-functions. 

Section 3 contains the crucial step in our analysis, showing that P n ,jv is an orthogonal polynomial 
with respect to an analytic weight supported on a contour in the complex plane. This leads us to a 
Ricmann-Hilbcrt formulation of P Ui n and the asymptotic analysis is performed in the subsequent 
sections via the steepest descent analysis [13] on the corresponding Ricmann-Hilbcrt problem. 

The proofs of Theorem 1.1, Theorem 1.2, and Theorem 1.3 are respectively in Section 4, Section 
5, and Section 6. 

Acknowledgements M. B. acknowledges a support by Natural Sciences and Engineering Re- 
search Council of Canada. S. Y. Lee acknowledges the support by Sherman Fairchild Foundation. 
K. D. T-R McLaughlin is supported in part by the National Science Foundation under grants 
DMS-0200749, and DMS-0800979. The authors wish to thank the faculty and staff of the Centre 
de Recherches Mathematiques for their support. The authors also thank AIM workshop "Vector 
equilibrium problems and their applications to random matrix models" and useful discussion in 
that research-conducing environment. 

2 2D equilibrium measure problem 

We introduce a compact region K as the support of the equilibrium measure p,* that minimizes 
C{p) defined in (1.4). From (1.5), it is well known that under certain regularity assumptions on 
the potential Q, the measure p* is known explicitly: dp* = ^ddQ ■ Ik- The central problem then 
becomes the determination of the support set K. 
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(2.1) 



For our specific potential Q given in (1.6) (which is of the type Q(z) = \z\ 2 — (harmonic)) 
^lfcdA is the unique equilibrium measure, if the following holds. 

(A) There exists £20 <= K such that U2B — on K; 

(B) For the £20 chosen above, U2D > on K c ; 

where 

Z&d(«) :=Q{z)- - / log\z - w\dA(w) + l 2D . (2.2) 

The following definition is specific 2 to Q given by (1.6). 

Definition 2.1. Given a compact set K C C 0/ area 7rt, (f t/iere exists a function S : K c —¥ C such 
that the following three conditions are satisfied 

i) S(z) has continuous boundary value on dK and satisfies S{z) = z on dK , 

ii) S(z) — c/(z — a) is analytic in K c , 

Hi) S(z) = {c + t)/z + 0(z- 2 ) for z -> oo, 

then we say that S is the Schwarz function of K or K has the Schwarz function. 

If K has the Schwarz function S, the following identities make a useful connection between the 
Schwarz function and U2B ■ 

{1 f wdw 1 f S(w)dw c 
z-\ * =z^ d> -^—^ = z , zGK, 
2i7r JdK z ~ w 27rl JdK z ~ w z ~ a ^ 3 ^ 

2i7r J q K z — w 2ni J q K z — w z — a 

Here the contour dK is positively oriented with respect to K. This identity gives 

dU 2 n(z) = {°J ZEK u (2.4) 

\z—S(z), z f K. 

We denote d — \{d x — \d y ) and d = \{d x + id y ); except when d is followed by a set (then it means 
"boundary"). From the equation (2.4) we obtain, for each connected component of K c , 



Hsd(z) = \z\ 2 - \z \ 2 - 2K 



S(w) dw 



z E (a component of K c ), (2-5) 



where zq G dK is any point on the boundary of the same component. 

Lemma 2.2. Suppose K is a compact set such that each connected component of K c is bounded by 
a finite union of closed real analytic curves (without any singularity) and K has a Schwarz function 
as in Definition 2.1. Then K is the strict local minima oflA2D i- n the sense that, for e > small 
enough (i.e. near the boundary dK), 

U 2 d(z + e(n x + in y )) = 2e 2 +0(e 3 ), Vz G dK , (2.6) 

where (n x ,n y ) is the unit vector that is outward normal to dK at zq. 



2 The standard notion of the Schwarz functions can be found in [12] 
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Proof. We use the fact that the Schwarz function S(z) can be analytically continued in a neighbor- 
hood of a real analytic curve dK. This then gives a smooth extension, say U (later at (2.24) we 
define this quantity as the maximal extension), 0HA2B in a neighborhood of dK by using (2.5). We 
have U = U 2 d on K c . 

On dK, dU(z) = dU{z) = 0. Because AU = 4, the Hessian has trace 4. Because |S"(z)l — 1 I0r 
z G dK (since S(z) = z there), the determinant of the Hessian of U, which is given by 

4 ([ddU{z)f - (d 2 U(z))(d 2 U(z))) = 4 (1 - |S'(z)| 2 ) , (2.7) 

vanishes. So the Hessian of U has eigenvalues and 4. Since the flat direction is tangent to dK the 
direction normal to dK gives the steepest ascent of U as described in the lemma. □ 

Lemma 2.3. Suppose K is a compact set of area irt such that each connected component of K c 
is bounded by a finite union of closed real analytic curves and K has a Schwarz function as in 
Definition 2.1. Then using this set K in (2.2), (A) in (2.1) is satisfied. If, furthermore, IA2D has 
no more than one critical point in K° then (B) is also satisfied. 

Proof. One shows (A) of (2.1) by direct calculation using Stokes' theorem and dU2n = dU2r> = in 
K by (2.4). 

Suppose there is a local minimum of IA2D at K c . Then the Mountain pass theorem says that 
there must be a critical point on a path connecting K and the local minimum. To briefly explain, 
one crosses a critical point z cr it (i.e. a mountain pass) by following a path for which the maximal 
elevation is minimal, i.e. U2T>(z C rit) == inf( pa ths)[ niax ze(path)^2D(-2 : )]-) This requires at least two 
critical points in K c . Therefore, if there is no more than one critical point in K c then there is no 
local minimum in K c , and the lemma follows. □ 

Plan of the proofs: In the rest of this section, we define the region K explicitly. Then we show 
that i) there exists the Schwarz function of K and ii) U2D has no more than one critical point in 
K c . By Lemma 2.3 this is enough to conclude that K is the global minima of U^d- 

Proposition 2.4. For t > t c , the global minimum ofU2D is exactly at the doubly connected region 
K = D(0, \Jt + c) \ D(a, t/c) where D(a, r) stands for the disk of radius r centered at a. 

Proof. Our proof relies on the Schwarz function of K: 



siz)= \(t + c)/z, \z\>Vt + c, 

\a + c/(z — a), \z — a\ < \fc. 

To show that there is no critical point of W2D at K c note that any point that satisfies S(z) — z is 
on dK. The proof is complete by Lemma 2.3. □ 

For the pre-critical case, the region K is defined by the following lemma. 

Lemma 2.5. For < t < t c , there exists a unique simply connected compact region K C C and a 
conformal mapping / : C \ D — > C\ K given by 

f(v)=pv peK+, cel + , a€(0,l), (2.9) 

v — a a 
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with its inverse conformal map denoted by F : ( £ \ K — >C\D, such that the function S{z) given by 
S(z) = f(l/F(z)), is the Schwarz function of K. The three parameters satisfy 



l-a 2 >-. (2.10) 
P 

The proof of this lemma is given in Appendix A. 

Note that the lemma asserts the existence of functions p(t), a(t), k(£); it is a direct consequence 
of the proof that these are analytic in t for < t < t c . 

In what follows we will need the deck transformation, defined via 

deck(u) := a- (2.11) 
p v — a 

This function is characterized as the unique nontrivial function satisfying f(v) = /(deck(u)). 

Proposition 2.6. For the pre-critical case, K defined by the conformal mapping in Lemma 2.5 
satisfies (A) and (B) of (2.1), i.e. the global minima ofhi2D are exactly at K. 

The existence of the Schwarz function follows from Lemma 2.5. We will complete the proof in 
Section 2.1, by showing that there is only one critical point of ^2d- 



2.1 Definitions of y(z) and F(z); Proof of Proposition 2.6 (pre-critical) 

There are two critical points of / : C — > C at a ± i\J nj p. By (2.10), both of them are inside D. Let 
us denote by vg the critical point in the upper half-disc, and let us denote by 8 the corresponding 
branch point in the z plane: 

8 := / I a + it /— I = ap h 2i^/Kp and vg := a + i, — . (2-12) 

V V p) a V p 

The branch point /3 has the positive imaginary part 2iy / Kp, and the other pair is given by the 
complex conjugates of the above. 

Lemma 2.7. The branch points (3 and f3 are in the interior of K. 

Proof. Since vg € D is a critical point of the 2-to-l branched covering map, vp is the only point that 
maps to (3 under / : C — > C. If 8 is not in the interior of K, then by the surjectivity of / : D c — > K c 
one of the preimages of 8 must be outside the unit disk, which is a contradiction. □ 

Because [v — a)(f(v) — 8) is a quadratic polynomial in v with double roots at vp we have 
(v — a)(f(v) — 8) = p(v — vg) 2 and, similarly, {v — a)(f(v) — 8) — p(v — vp) 2 . We also have 
(v — a)(v — deck(w)) = (v — vp){v — vp) which, combined with the above, gives 

(t> - deck(v)) 2 = p- 2 (f(v) - 8)(f(v) - J). (2.13) 

We also obtain by a direct calculation that 

1 , I /' 2 '•' x - 1 



v + dcck(v) = -f(v) + - [a z + - = -(/(«) + |/3|), (2.14) 
P a V PJ P 
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where we use 



(2.15) 



that is obtained by taking v — in (2.13). From (2.13) and (2.14), for any v ^ {vp, v@}, there is 
exactly one choice of sign such that 



1 

2p [ 



f(v) + \[3\±J(f(v)-p)(f(v)-f3) 



(2.16) 



is satisfied; The other sign gives deck(u). Now let us select our branch cut to be inside K and to 
intersect K only at the negative real axis. This is possible because K n R_ ^ which is simply 
confirmed by /(— 1) = — p — n/(a(l + a)) < 0. Having determined the branch cut, we select the 
branch of the square root by [(z — (3)(z — P)] 1 ^ 2 ~ z + 0(1) as z — > oo. Then (2.16) gives the 
following lemma. 

Lemma 2.8. Take any simple smooth curve B in K that starts at j3 and ends at j3 while intersecting 
R only once at the negative real axis. The map F (the inverse of f; see Lemma 2.5) has the analytic 
continuation over C with the branch cut at B via the formula 



F(z) 



1 

2f> 



z+\p\ + J(z-p)(z-p) 



, deck o F(z) = 



1 

2f> 



z + \P\-<J(z-0)(z-I3) 



(2.17) 



Proof. It is enough to determine the sign in front of ... by considering F(z) w z/p as z — > oo. □ 

Using this F, the analytic continuation of the Schwarz function is given by S(z) — f(l/F(z)), 
of which the explicit form is announced in the next lemma. 

Lemma 2.9 (Definition of y{z)). S(z) has analytic continuation to the entire complex z plane with 
the branch cut at B (that is defined in Lemma 2.8) via the formula 



S{z) 



a c 1 c + t 1 , . , , ( z - b )\J - P)i z - P) 
2 + 2—a + ^z--2 y(z) > V{z) ^ a (z^afz ' 



where y(z) satisfies three conditions respectively near a, and oo: 

— c 

y(z) , y(z 



— , y(z)~a--+0(l/z 2 ). 

z z 



(2.18) 



(2.19) 



We also note that the middle condition in (2.19) and the last condition in (A. 2) gives 

b=\p\/(t + c)=a/p. (2.20) 
It is useful to define a function S^ack (^ ne analytic continuation of S through the cut B) by 

SbackW := /(l/deck(F(z))) = S(z) + y(z). (2.21) 
The following Lemma, combined with Lemma 2.3, proves Proposition 2.6. 
Lemma 2.10. The only point that satisfies S(z) = z and F(z) ^ dD is at z = b G (a, oo). 
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Proof. If z satisfies z = S(z) then v = F(z) satisfies f(v) — f(l/v) = f(l/v). This implies 
deck(w) = 1/v if v ^ (3D, and we have 



S(z) = /(1/v) = /(«) = f{v) = f(l/(l/v)) = /(l/deck(«)) = S hack (z). (2.22) 

Therefore, z — b because Vp is fixed under the deck transformation. 

To show that b € (a, oo), one notes that KPiR. is connected (otherwise K is not simply-connected) 
and lies to the left of a (which can be checked from the conformal map). Because W2D grows to 00 
at a £ K and +00, U2D ■ K — > K must have at least one local minimum at (a, 00), which must be a 
critical point of U2D, where S(z) — z holds. □ 

The last lemma tells that $2_d has at most one critical point in K c , which proves Proposition 2.6. 
2.2 Branch cut of y(z): pre-critical case 

Having defined S(z), y(z), F(z), as the functions on the whole complex plane with the branch cut 
on B, we can define the following objects 



U OP (z) :=5R 



y(w) dw 







U(z) := |z| 2 -|z | 2 -25R 



S(w) dw 



(2.23) 



for a point z s.t. F(z ) G <9B, (2.24) 



where the integration contour is in C \ B. From (2.5) we note that U is the smooth extension of 
U 2 b from C\K. 

In the previous subsection, we have not specified the exact location of the branch cut B except 
that it is inside K and intersects the negative real axis. In this section we determine the precise 
location of the branch cut in such a way that will be useful to the subsequent Riemann-Hilbert 
analysis of the orthogonal polynomial. 

Definition 2.11. The level line(s) from a point, say p, is defined to be the connected component 
of the set {z : 5ft [/ y(w)dw] — 0} which contains p, and, since the integral defines a continuous 
function up to sign, the level line is given independently of the branch cut. 

We first prove that the level line from j3 does not intersect the level line from b. 

Lemma 2.12. We have 

2U OP (z) - U{z) +U(S back (z)) - \S back (z) - z\ 2 . (2.25) 

The proof uses the anti-holomorphic involution structure. We leave the proof in Appendix C. 
We obtain the following corollary that tells that the level lines from (3 and b do not intersect. 

Corollary 2.13. We have 5ft[ J„ y{z) dz] = Wg£i(6) > 0. The contour of integration is in C \ B. 

Proof. By Lemma 2.10, we immediately obtain the first equality. The last inequality is from Propo- 
sition 2.6. Note that, by Lemma 2.10, b is not in K so 14(b) = L(2v(b). □ 
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Figure 4: Possible topological configurations of level lines from (3 and b. The last (boxed) configu- 
ration satisfies the inequality in Lemma 2.13. The thick lines are the (potential) branch cuts. 



In Figure 4 we show all the three possibilities that the level lines from /? and b can place 
themselves in topologically distinguished ways. We explain this below. 
The rules for the level lines from b are as follows. 

{The lines are symmetric with respect to the real axis; 
A closed loop of level line must include pole(s) at and/or a; ^ 
Four level lines emanate from b; 
Out of four, at most two level lines can go to infinity; 

Following these rules, we first draw all the possible level lines from b (the dotted lines). If none of 
the level from b escapes to oo then the level lines from b follows the dotted lines in the left figure; 
If one pair of the level lines escape to oo then there are two possibilities, the middle and the right 
figures. 

With the dotted lines being drawn, the level lines from j3 are now uniquely determined. The 
rules for them are exactly same except the last two items must change from "four" to "three" . 



Out of the three possibility, we prove that only the last one satisfies the inequality in Lemma 




Figure 5: y(x) on R: schematic view. 



1G 



2.13. Note that y{+oo) = a > 0. Therefore y(x) > if x > b and y(x) < if a < x < b, 
because there is no branch cut intersecting on [0, +oo). (See Figure 5 for the behavior of y(x) on 
R.) Therefore, when we restrict x G (a, +oo), J, y(s) ds reaches the minimum at b. This means 
that, for the left and the middle cases in Figure 4, the dotted lines are "lower" than the regular 
lines and 5ft [/a y(z) dz] < 0, which contradicts Lemma 2.13. Therefore we have proven the following 
lemma. 

Lemma 2.14. For all < t < t c , the level lines from b as well as those from (3 and (3 are as 
shown in the last (boxed) configuration of Figure 4- More precisely, the level lines from b consist of 
three analytic arcs, one which encircles the point a and two which diverge to oo. Two components 
of the level lines from f3 are analytic arcs which terminates in (3, and together form a closed curve 
encircling 0. 

We want to choose the branch cut of y(z) from one of the level lines from (3. For this to be 
consistent with the previous definition of B (that appeared in Lemma 2.8), we only need that the 
chosen level line be inside K. Below we prove exactly this. 

The logic requires to use lots of hats (""") temporarily. We define B by the unique level line from 
(3 to f3 that intersects the negative real axis. Then we define the hatted version of the functions: 
S(z), y(z), F(z), by the formula (2.17) and (2.18) with the branch cut at B. After showing that B 
is inside K, we can replace all the hatted objects by the un-hatted ones. 

Definition 2.15. We define B by the unique level line from (3 to (3 that intersects the negative real 
axis. We place the branch cut of S(z) and y(z) exactly at B using the formula (2.17) and (2.18). 
We also use (2.23) and (2.24) accordingly to define Uop andlA? 

Lemma 2.16. For any p G B\{(3, (3} let us define the normal derivative 9 n := n x d x +n y d y atp along 
either of the two unit vectors n := (n x ,n y ) that is perpendicular to B atp. Choose + and — sides of 
the cut such that n points toward the + side. Then we have d n ltop{p)\ + = — d n Uop(p)\_ = — \y(j>)\- 

Proof. The absolute value of d n Uop(p)\ + is \y(p)\ because d n ttop(p)\ + is purely real. Defining 
n = n x + in y we have (the other — side works similarly) 



d n U p(p)\ + = (nd + nd)Uop(p)\ + ^ s H[n-y(p)\ + ) . (2.27) 

By the definition of B, n ■ y(p) is purely real and either of ±|y(p)|. The correct (— ) sign can be 
determined at the point B n R from Figure 5. □ 

Lemma 2.17. There is no local minima ofU{z) on B. 

Proof. Taking the same notations as in Lemma 2.16 we have, for p G B\ {(3, /?}, 

[dnU(p)] + - [£>„%)] _ = [dMopip)] + - [dMovip)] _ = ~2\m\- (2-28) 

This means that at least one of [<9 n Z//(z)] + and — [<9 n W(z)]_ is negative. Therefore, U(z) decreases 
in one of the directions perpendicular to B. 



3 The existence of a point zo follows from {j3,/3} C K: starting from B that is completely inside K (so that all 
z £ dK satisfies S(z) = z) one cannot make B cross the whole of dK by any continuous deformation B — > B. 
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If p = ft we consider the expansion around ft: 

U(z) = (real analytic in R 2 ) + Uop(z) 

= U(ft) + Cl (x- 3?/3) + c 2 (y - 3/3) + 5R(c (z - ftf' 2 ) +OQz- /3\ 2 ) ' ' 

where c\,c% 6 K and cq € C is determined to match y(z) = fcoV 1 z — ft + 0({z — ft) 3 ^ 2 ) and Co 
is nonzero. If either of C\ or c-i are non-vanishing, then there exists a point z near ft such that 
< Li (ft). If both ci and C2 vanish, there still exists such a point that lowers the value of U by 
the non-vanishing term $1(cq(z — /3) 3 / 2 ). This proves Lemma 2.17. □ 

Lemma 2.18. For the pre-critical case, we have B C K , and, by taking B — B, U(z) has the global 
minima exactly at dK . 

Proof. If B intersects dK at p then any neighborhood of p contains a region where IA = U, and 
U(p) = U(p) by continuity (U is continuous by definition). By the previous lemma, p cannot be the 
global minimum of IA and there must be a point z ^ B where U reaches the global minimum. This 
point satisfies \F(z)\ ^ 1 because, for \F(z)\ — 1, one can show that U{z) = U(p). Then z must 
satisfy dli(z) = z — S(z) = and, by Lemma 2.10, z must be b. Both B and B do not intersect 
[0, +oo) and we have U[b) = U{b), and U{b) -U{p) = U(b)-U{p) > by Corollary 2.13. Therefore 
B does not intersect dK, and dK is the global minima of U. □ 



2.3 Branch cut B: The post-critical and the critical case 
Post-critical case: Wc define y by 

/ c t + c\ (z-b)(z-f3) (+ : z e Ext(B), , n on , 

y(z):=± U+ )=±a± M; , <^ (2.30) 

V z — a 2 / — a) I - : z e Int(B), 

for some closed Jordan curve 6 that we determine soon. In fact, this formula is consistent with the 
one (2.18) for the pre-critical case, if one uses the two formula in (2.8) respectively for Ext(2?) and 
Int(S). Both ft and b are uniquely determined by the equation above as (they are both real and we 
make a choice of ft < b) 

u _ a 2 + t+^(t-a^-^c a _ a 2 +t-y/(t-g2)2-4a2c 



One can check that a is less than ft. With such configuration, there is topologically only one 
way to draw the level lines from ft, which is illustrated in Figure 6. A proof is similar but simpler 
than that for the pre-critical case. Consider all the (three) possible ways of drawing level lines from 
b using the rules in (2.26); there exists only one case where the level lines from ft can be drawn 
consistently with the level lines from b. 

Critical case: For t-t c = 0{N~ 2 ^) we choose y(z) to be given by the equation (2.30) except 
that we extend the definition for t < t c . As a result, the two critical points y(z), b and ft, are exactly 
given by (2.31) (taking the principal branch of the roots). Their asymptotic behavior is given by 

c l/4 c l/4 

b = b c + — r s/t=t c + 0{t-t c ), ft = b c - —=Vt^T c + 0(t-t c ), (2.32) 
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Figure 6: Topological configuration of level lines from (3 and b for the post-critical case. 



with b c := a + y/c. For t > t c the level lines from j3 are exactly as in the post-critical case. For 
t < t c , one can verify, similarly as in the post-critical case, that the level lines from /? is as shown 
at the left side in Figure 8. We skipped the proof. 

Definition 2.19 (B, y(z), U, Uop)- For the critical and the post-critical case, we define B by the 
unique level lines of (3 that is a Jordan curve enclosing both and a. We define Ext(£>) and Int(£>) 
respectively to be the exterior (including oo) and the interior of the Jordan curve B. We define y(z) 
by (2.30) and then S(z) is accordingly defined over the extended domain C \ B by 

Using the above defined y(z) and S(z), we define Uop andlA by (2.23) and by (2.24) (where we can 
choose explicitly zq = y/t+~c). 

We note the reader that, for the post-critical case, B is not inside K, and S(z) does not match 
S(z) at (2.8) in some part of {z : \z — a\ < y/c}. 

Lemma 2.20. For the post-critical case, B sits in P(K) = {z : \z\ < \fi + c}. Also lA(z) has the 
global minima exactly at &P(K) = {z : \z\ = \Jt + c}. 

Proof. Since the proof is similar to the pre-critical case, let us skip some details. One can show 
that the global minima of U cannot occur onB \ {/3} as in the proof of Lemma 2.17. One can show 
(i G B is in {z : \z — a\ < y/c} by using yit — a 2 ) 2 — Aa 2 c > t — a 2 — 2a\[c with the formula for 
(3 at (2.31). The point j3 is ruled out from the global minimum by explicitly comparing U(f3) with 
U(\/t + c). Therefore, the global minima can only occur in {z : \z\ = yr + c} U {z : \z — a\ = ^Jc\ 
(those are the only places where dU(z) —z — S{z) — can possibly occur), and B \ {/?} should not 
intersect both circles at the same time. This means that B is in the interior of the larger circle and 
(whether B intersects the smaller circle or not) the global minima of U is on {z : \z\ = \/T+~c}. □ 

Lemma 2.21. For t = t c , B \ {(3} sits in K = {z : \z\ < y/t + c and \z — a\ > y/c}. Also U(z) has 
the global minima exactly at dK = {z : \z\ = \fi + c} U {z : \z — a\ — \/c}. 

Proof. By the same logic as in the pre-critical case, there is no global minima on B\ {/?}. Then the 
only possible global minima is where S(z) = z, which is the big circle and/or the small circle. Both 
are indeed the global minima because they are the limit of dK as t /* t c where dK is the global 
minima for all t < t c by Lemma 2.18. As a consequence, the cut B \ {(3} is strictly inside K. □ 
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3 Orthogonal polynomial on a contour 

As the first step toward the asymptotic analysis of P n .N we show that P n ,N is an orthogonal 
polynomial with respect to a non-hermitian, n-dependent orthogonality on certain contours in the 
complex plane. This observation hinges on the following lemma which relates the area integral to 
a contour integral. 

Lemma 3.1. The following identity holds 

JJ j T^a) k \z - a\ 2Nc e- N ^dA(z) = I z j w k+1 , N (*) dz, (3.1) 

c 1 r 

where 



(z - a) Nc e~ Naz 

™ n ,N(z)-=- ^iWn ' ^ C \ [0, a] . (3.2) 

The contour of integration T is a closed loop enclosing both and a counterclockwise. 
Proof. The integrand in the left hand side of (3.1) can be re-expressed using the identity: 

dx(z) = {Y^) k \z - a\ 2Nc e- N ^\ X (z) ■= {z ~ a) Nc [* (s - a) Nc+k e - Nzs ds (3.3) 



where z G C \ (— oo, a) and the contour of integration is chosen so as not to intersect the half-line 
(—oo, a). 

We claim that the continuous extension of x(z) is smooth in C. Let us show that the boundary 
values on both sides of (— oo, a) match. Let us take x < a and use a notation x± — x ± iO. The 
boundary values of the integral on both sides are related by 

" (s - a) Nc+k e- Nxs ds = e- 2 -( Wc + fe ) f ~ (s - a) Nc+k e- Nxs ds (3.4) 



Due to the prefactor (z — a) Nc , x( z ) 1S single- valued and continuous across (— oo, a]. 
The integral in x( z ) can be divided into two terms by 

r l>~ZXOO p~zxoo 

(s-a) Nc+k e- Nzs ds= (s-a) Nc+k e- Nzs ds- (s - a )^+k -Nzs r 

T(Nc + k + l)e- Naz 



{Nz 



0(e- N \ z \ 2 ' 2 ) , 



(3.5) 



where the branch cut of z Nc+h+1 lies on (— oo,0]. 

We derive the last equality below. For fixed N, c, k, there exist R > such that the above 
integrand satisfies Us — a) Nc+k e~ Nzs | < \ e Nzs / 2 e ~ Nzs | for all s € \z, z x oo) and \z\ > R. Therefore 
the last term in the first line in (3.5) is bounded by 



( s _ a )» c+k e-" zs ds 



< e- N "^ds = 2 , \z\>R. (3.6) 

Jr n \z\ 
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The first term after the first equality in (3.5) becomes, by the change of variables: X — z(s — a), 

rzxca -Nza f+oo p/ N , , , -, \ Naz 

/ (s - a) Nc + k e- N * s ds = 4™ / X Nc+k e- NX dX = 1 + * + . (3.7) 

J a V ' z Nc+k+l J q (Nz) Nc+k + 1 x ' 

Since \ is smooth, we can apply Stokes' theorem to get 

(LHS of (3.1)) = lim f f z^z-af \z - a\ 2Nc e~ Nz *dA(z) (3.8) 

r— too J J 
\z\<r 

= lim ff z 3 dx(z)dA(z) = lim — <£ z J x{ z )dz (3-9) 

11 ^(■-< f |y ffr° ^(e-^'')) fc P-10) 



lim 



The first term is independent of r as long as the circle \z\ = r encloses and a; and the last term 
vanishes in the limit r — > oo. □ 

This lemma allows us to view P Ui n as an ordinary orthogonal polynomial: indeed the system 
(1.1) is equivalent to 

( 7 _ n \Nc-Naz (0 , k = 0, 1, ••• ,n- 1, 

^.A^ r J e+fc+1 dz = jV^+»+l (3.11) 

l 2 V(iv c + n + i) ^' fc = n - 

By relabeling fc — > n — 1 — j, we get a standard system of non-hermitian orthogonality relations 

. f z _ a )Nc -Naz 

P n,N(z)z^ -^- n dz = (j = 0,l,...,ra-l). (3.12) 

Remark. The proof of Lemma 3.1 can be extended to treat weights of the form 

K 



W(z,z)=e- N ^ 2+?li ( az2+ ^l[\z 



aA 



and allows to relate the bimoments J c z J z k W(z,z)d 2 z with the moments of an iterated contour 
integral; this realizes the orthogonal polynomials as the biorthogonal polynomials appearing in [4] 
for a specific semiclassical bilinear moment functional. In particular one could use the formalism 
in loc. cit. to produce a corresponding Riemann-Hilbert problem (whose dimension depend on the 
integer K, however). We do not pursue this issue here. 

3.1 The Riemann-Hilbert Problem 

It is well-known that the polynomial P Ui n(z) can be characterized by a 2 x 2 Riemann-Hilbert 
problem [19] which we recall in our setting. Let us define 

Vk ■= j> z k w n , N (z)dz. (3.13) 
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Introduce the polynomial 



i-l,JV 



dct I 



■ det 



I +jJ0<ij"<n-l 



Vl 



v n -i ■ ■ ■ 
1 z 



v n -i 



u 2n-3 



(3.14) 



Note that Q n -i,N is not necessarily monic and its degree may be less than n — 1: its existence relies 
uniquely on the nonvanishing of the determinant in the denominator, which will be shown below in 
Proposition 3.2. It is also orthogonal to the powers 1, z, . . . z n ~ 2 . Define the matrix 



Y{z) = 



Pu,n(z) 
-2mQ n _x N {z) 



1 

2~7ri 



Pu,n(z') 



Q, 



-l.N 



CO 



/r z' — z 

Then Y satisfies the standard [19] Riemann-Hilbert problem: 

Y(z) is holomorphic in C \ T; 

1 Wn iN {z) 



w n< x{z')dz' 
w n , N (z')dz' 



(3.15) 



Y + (z) = Y_(z) 
Y(z) =(l + 



zeT: 



(3.16) 



if and only if det [fi+j] <j j< n -i ^ ; this is seen by noticing that if Q n ,N exists as given by (3.14) 
then 

Qn - 1 ^ z '\ WntN { z ')dz' ~ z~ n (l + Oiz- 1 )) (3.17) 



it z' — z 

We thus need to show the nonvanishing of the determinant in the following proposition. 
Proposition 3.2. The determinant det^i+^JcKijXn-i does not vanish. 
The proof can be found in Appendix D. 

3.2 The g-function and the 1st transform 

Given the Riemann-Hilbert characterization in eqs. (3.16), this asymptotic case can be addressed 
using the Deift-Zhou [13] steepest descent method as soon as an appropriate g-function has been 
constructed. 

The problem here is that we cannot rely upon the minimization of an energy functional as in the 
case of ordinary orthogonal polynomials on the real line. Instead, we will define g-function from 
y(z) that is defined at (2.18)(2.30) and in Definition 2.19 for all three cases. 

The weight function can be expressed in terms of a potential V (of the semiclassical type [23]) 

by 

w n , N (z) = z~ r e~ NV( - z \ V(z) :=az~cln(z-a) + {c + t)\nz. (3.18) 
We choose the cut of the logarithms at [0, +oo). 
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Lemma 3.3 (Definition of Tbp and rW). The steepest "ascent" path (the maximal elevation) of 
Uop from j3 reaches b. Equivalently, the steepest descent paths from b reach f3 and /3; Let us denote 
the former (that reaches f3) by Ybf) and the latter by r^. 

Proof. In the last part of Appendix C we show that J~ y(s) ds is purely real. This means that 6 is 
on the steepest ascent path from (3. (See Figure 7 for a numerical confirmation.) □ 




- 1 1 2 



Figure 7: B (thick line); level lines (dotted lines); steepest descent/ascent paths (the other lines). 
All the lines are with respect to the function Uop and emanate from (3 or b. We close the contour T_ 
at some region 5R(z) < (by deviating from the steepest descent paths) such that £!_ is a bounded 
region. (All the lines are numerically plotted for t = c = a = 1.) 



Definition 3.4. For all t, we define <f>{z) by 

y(s)ds (3.19) 



where the integration contour lies inC\(6U T^p U \b, +oo)) for the pre-critical case and raC\(6U 
[/3,+oo)) for the other cases. We haveU p{z) = 3?0(z) by (2.23). 

Definition 3.5. For all t, the g-function is defined by 

2g{z):=jV{z)-- t <t>{z)+L (3.20) 
The constant £ S K is defined such that g(z) — log z — > as z — > oo. 
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Figure 8: The level curves and contours for the critical case with t < t c on the left; and t > t c (for 
both the critical and the post-critical case) on the right. The branch cuts B (thick lines) are chosen 
by one of the level lines from /}; the other level lines from j3 and b are shown with dotted lines; the 
other lines are steepest descent/ascent paths. (All the lines are numerically plotted for c = a = 1 
and t = t c - 0.02 (left) and t = t c + 0.02 (right) where t c = 3.) 



Lemma 3.6. For the pre- and the post- critical cases (not for the critical case), we have 

1 



9{z) = 2nt J B l0 S( z ~ s )\y( s )W ds \ ( 3 - 21 ) 

where we take the branch cut of the logarithm in B U T^p U [6, +oo) for the pre-critical case and 
BU [/3, +oo ) for the post-critical case. 

Proof. One can check that g'(z) is analytic away from B; the singularities of y(z) as described at 
(2.19) exactly cancel those from V'(z). Let us only consider the pre-critical case as the post-critical 
case is similar. To check the jump discontinuity of g on z € B we evaluate 



1 

2~t 



z + 



y(s) ds 



1 
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y(s)- ds 







/ y(s)+ ds 


-il 


z + 



\y(s)\\ds\ (3.22) 



where ± sides are with respect to the orientation of B directed from /3 to /3 for the pre-critical case 
(and counter-clockwise for the post-critical case). The integration contour is taken away from the 
cuts of both V and <f>. In the last equality, we use that y(s) ds = ±z|y(s)| |ds| for the infinitesimal 
line segment ds tangential to B. The sign can be determined at the intersection of B with the real 
axis using Figure 5; as in the proof of Lemma 2.16. The above calculation gives exactly the jump 
of the proposed formula for the ^-function. □ 

Lemma 3.7. On the exterior ofP(K), the g-function is the complex logarithmic potential from 
the probability measure that is uniformly supported on K, i.e. 

1 



Ug(z) =3? 



log \z — w\ dA(w) 



K 



z £ P(K). 



(3.23) 
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Above, the exponential is applied to disregard the (unimportant and arbitrary) branch cut. 
Proof. It is enough to show that 



g'( Z ) = L(y'(z)-y(z)) = ^[S(:,- 



= l 9 {Q{z)-U{z)) = - 

t TT 



dA(w) 



K 



(3.24) 
□ 



where each equality is from (3. 20), (2. 18), (2. 24), and (2.2), in the order of appearances. 

From (the proof of) Lemma 3.6 the sign of dicf> is negative on both sides of the cut B. 
We may choose the contour T (of orthogonality) by any continuous deformation within C\ [0, a 
We choose 

B U r &( 3 U Y b -p pre-critical, 
B critical and post-critical. 

We write some useful relations that follow from Uqy> — 5i</>. 



r := 



(3.25) 



l/z r , 

„-NUop 



zeB; 
zeT\B. 



(3.26) 



For the critical and the post-critical cases, the ^-function defined at (3.20) turns out to be 
simple: 




z e Ext(S); 
z e Int(B); 

with Robin's constant £ is given by (such that to make Sfti? continuous at (3) 

I = SR [((< + c ) log P - c log(/3 - a) - /3a) ft] . 



(3.27) 



(3.28) 



Definition 3.8 (T± and f2±). We choose contours T + and T_ following the steepest descent paths 
(as shown in Figure 7 and Figure 8) from (3 (and j3 for the pre-critical case) such that 5R0(z) = 
Uop(z) < on those contours. We choose T_ to be the unbounded ones. We then deform T_ 7 away 
from the steepest path, to make it a single bounded curve such that IAop{z) < on the deformed 
curve. Then the domains fl± are defined by the open sets enclosed by B and T± respectively. See 
Figure 7 for the pre-critical case and the right side in Figure 8 for the post-critical case. 



1 w n . N (z) 
1 



On T (3.25), we have 

Y + (z) = Y_(z) 
Correspondingly we introduce the matrix A{z) as follows: 

/ zec\o + un_ 



(3.29) 



A{z) := e 



= P. 2 



'Y(z) < 



1 

-l/w ntN 1 
1 " 



z e fL 



(3.30) 
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The RHP for A is 



A + = A. 



A,=A_ 



A+ = A- 



1 " 








z~ r ' 




-z r 




z -r e -NUo? - 


1 



A(z) = [I + 



z 




on T±; 



on B: 



on T \ S; 



z — > oo. 



(3.31) 



The remaining jump of g is ±2in and hence e~ tN 9i z ) does not have additional jumps besides the 
ones indicated above. On T±, Uop = $t4> is negative, and on T \ B (pre-critical only) it is positive 
and uniformly so except near the endpoints /3 and /3. At any finite distance away from 6, the jump 
matrices in (3.31) are exponentially close to the identity jump in any LP norms and in L°° norm. 
As the nonlinear steepest descent method goes, one first constructs a model RHP where only the 
jump exactly on B is retained. This is done below. 



4 The pre-critical case: t <t c 

Outer parametrix: Wc define the following model RHP for <I> by 



§4 



1 + 





-z r 

oi 1 



on B , 








(4.1) 



There are several ways of expressing the solution $ and we provide the one that is most relevant 
to the geometry of K. 

We use / : C — > C as the uniformization mapping of the domain 
with the branch cut C \ B. This map is 2-to-l branched covering 
map and we have defined the inverse map F at (2.17). There are 
regions, E = F(C \ B) and E = deck o F(C \ B), such that E and E 
are each conformally equivalent to C \ B. We define j± = [F(B)]±, 
the two preimages of B under / when approached from ± sides. See 
Figure 9. 



Proposition 4.1. The scalar function given by 



D(z) 



F(z) 



(4.2) 









7-T 


s 


7+ 











satisfies the following properties: 

1. D(z) is analytic and nonzero in C \ B 



Figure 9: Uniformizing plane; 
See the text. 
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2. D + (z)D-(z) = z, ze B; 

3. lim z _ > . 00 D(z)/ z > 0. (Let us call this positive number by D^.) 
Proof. The proof relies on the identity 

f(v) = (p / a) v deck(v) . 

This can be shown by 

= (v — a)(deck(u) — a) + n/p = v dcck(w) — a(v + dcck(v)) + a 2 + n/p 
= v deck(w) — (a/p) f(v), 



(4.3) 



(4.4) 



where the last equality comes from v + deck(v) — (f(v) + \/3\)/p at (2.14) and |/3| = (p/a)(a 2 + n/ p) 
that is shown right below (2.14). Using (4.3) and that F{z)\ + = deck(F(z))|_ on z € B, we have 



D + {z)D(z) = - deck(F(z)) • F(z) = f(F(z)) = z, z e B. 
a 

The last condition is checked by — lim z _ >00 \J p/a F(z)/z = 1/^/pa > 0. 
Proposition 4.2. The solution to (4-1) is given by 

{paYinpY' 2 



(4.5) 
□ 



$(z) = y/ P F>(z) 



1 

<pa)- r ^p) 1/2 
pF(z) — pa 



pF(z) - pa 
1 



(pF(z)Y 



(4.6) 



We choose the principal branch of the square root, i.e. \J pF'(z) = 1 + 0(1/ z) as z — > oo. 

Proof. Though the solution may be checked by a direct calculation, here we present a self-contained 
derivation (which is essentially the spinor construction in [6]). 

For the solution there exists an analytic matrix valued function H:S-> SL(2,C) such that 



Il(F(z)) = (D oo y^^(z)D(z)-^. 
This satisfies the RHP (4.1) of $ for r = 0, i.e. 



U(F(z))+=U(F(z)). 



1 

-1 



H(F(z))=I + 0(l/z), 



We may freely extend H on the whole complex plane by putting 



H(«) := H(deck(t;)) 



1 

-1 



(4.7) 



(4.8) 



(4.9) 



This makes H to be continuous on 7 + , by the jump condition in (4.8). On v e 7_ we evaluate the 
jump as we approach from S and E by 



n(v)L = H(dcck(v)) 



-1 

1 



-H(u)|~, v G 7- 



(4.10) 
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The last equality uses (4.9) with the continuity of H(deck(w)) across 7_. Then RHP for H is given 

by 



H(v)L = -H(u) ? 



H(«) = 1 + 
H(«) = 



1 

w 

1 

-1 



(4.11) 



+ O (v — a) v — > a . 



This gives a scalar RHP for each entry; One may only solve the (ll)-entry and the (22)-entry as 
the others are then given by (4.9). We obtain 



H(«) 



4/2 



i 



(wp) 1 / 2 v — a 
where we take the branch cut on 7_ such that 

,1/2 



1 K 
(k/o) 1 / 2 u — a 
1 



V/'M = — J (v ~ vp)(v - vp) « <j P ' 

u — a V I k l > /(v — a), v 



(4.12) 



u — > oo; 

Q . 



(4.13) 



Using f'(F(z)) F'(z) = 1 and the formula for L>(z) at (4.2) we obtain to be 

1 

1 — 

*(z) = (pa)5 CT Vp^'(^) 



1 



(ftp) 1 / 2 F(«)-a 



(rep) 1 / 2 F(z) - a 
1 



F(z) 



This proves the proposition. 



(4.14) 



□ 



Local Parametrices: The procedure for constructing the local parametrices near the points j3, (3 
is standard and the reader may refer to [13, 14]. Here we only consider a neighborhood around /3 
and skip the similar construction around j3. 

Let D» be a small disk centered at j3 with a finite radius. We define the local coordinate C by 



C(z) 3/2 := N (4>{z) - m) . z G 



(4.15) 



such that £ maps P 6 -g fl D-g- into K+. Then £ maps T± respectively into the rays [0, e ' oo). 
This defines a conformal map £(z) from the fixed disk to a domain in the £-plane that expands as 
0(N 2 / 3 ). 

We want to find the function Vj(z) in (and, similarly, Vpiz) in D^) such that 



{satisfies the jump of A at (3.31); 
converges to &(z) on the boundary of as N — > oo. 



(4.16) 
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This is satisfied by the following Ricmann-Hilbcrt problem: 



-1 

1 



^(z)~/ + 0(C(z)- 3/2 ) 
V-Az) ~ G(C(z)- 1/4 ) 



1 

-1 



1 

elCC*)"" i 

1 e -S«*> 8/1 
1 



on r± n D s 



on r b ^ n 



as z 



as z — > j3 



(4.17) 
(4.18) 

(4.19) 

(4.20) 
(4.21) 



The last condition is necessary later for "the error matrix" to be analytic at f3. 

The solution to the above RHP is given in terms of the standard Airy parametrix A(() by 



e 4 



1 

71 



i -i 
i i 



where 



AO ■- 



2ne 4 < 



2/o (0 "2/2 (C) 

y'o(C) -y' 2 (0 

-vi (0 -2/2(0 
-2/U0 -y 2 (C) 

-2/ 2 (C) 2/1 (0 
-2/ 2 (0 2/i(0 



2/o (0 2/i (0 
2/o(0 2/UO 



C(*)^(C(*)) , 

K*os argCG (0,2tt/3); 
e§£ S<T3 argCe (27r/3,7r); 
fC5<r 3 argC e (7r,57r/3); 

argC G (5vr/3,27r). 



(4.22) 



(4.23) 



We have used the standard Airy function to define: 

yj (C) :=^'Ai(^'0, J =0,1, 2, w 



- „2i7r/3 



(4.24) 



Using $ (4.6) and Vj (4.22) we define ^4°° which we show (in the following Error analysis) to 
be the leading approximation of the matrix A(z) defined at (3.30). 



!$(z), for z outside of the disks Bp WDp 

$(z)z~( r / 2 ) <T3 7> 3 (z)z( r / 2 ) <T3 , inside D^; 

$(z)z-< r / 2 )' T3 ^(z)z( r / 2 )^ ; inside %. 

Error analysis: We define the error matrix £ by 

£{z) := A°°{z)A- 1 {z). 



(4.25) 



(4.26) 



As a consequence of the jumps of A and A°°, £ has jumps only on the boundaries of the disks, 
and D» and on the parts of T± and T\B that lie outside of those disks. 
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The jumps of £ on the boundary of the lens, r± \ (I5g U D^), is given by 

[£{z)] + [{£)-\z)]. = [A°°{z)] + [A-\z)] + [A{z)]-[{A°°)- l {z)]- = *{z) {l + 0{e^ N )) *-\z) 

= I + 0{e- cN ) , for some c> 0. (4.27) 

The jumps of £ on the boundary of the disc, dH-g, is given by (plus side is inside the disc) 

[£ (z)}+[(£)-\z)}- = *(z)z-WV**V^{z)z { rM' r '*-\z) = 1 + ©(A^ 1 ) , (4.28) 

using the boundary behavior of (4.20) and (4.15). So £ satisfies a small-norm RHP and we conclude 
that 

' 1 



A(z) 



I + O 



N 



A°°(z) 



(4.29) 



uniformly in (a compact set of) the (extended) complex plane. 



Strong asymptotics: The asymptotic for the orthogonal polynomial can now be read off from 
our approximation (4.25) by tracing back the transformations from Y{z) (3.15) to A(z) (3.30). 
Though the computation is straightforward we show a few steps for the interested readers. 



Away from B: The asymptotic for the orthogonal polynomial can now be read off from our 
approximation using that P u ,n{z) — \Y{z)]n — A\ie tNa ^ outside of ft±. Using the asymptotic 
behavior (4.29), we get the strong asymptotics in C \ (D^ U O-g- U f2 + U f2_) as follows. 



P n , N (z) = [(I + 0(1/N)) A°°] u e tN ^ 

= (1 + 0{1/N)) ^F^)(pF(z)Ye tN ^ 



(4.30) 
(4.31) 



Near B but away from (3 and j3: Tracing back the transformations from Y(z) (3.15) to A(z) 
(3.30) we find, respectively for z e fi± \ (D^ U Dj), 



Pn(z) 



e tiV|<T 3 



A(z) e 



A(z) e 



tN(g(z)-^)a 3 



1 

±l/w n N 1 



11 





1 " 


e f V(z)a 3 




±z r 1 





(4.32) 
(4.33) 



Recalling that A(z) = (I + 0(N 1 ))A co {z) and A°° = $ in the region of our interest, we find 
respectively for z £ f2± \ (D^ U !>§•), 



P n (z) 



- „f(v(z)+a) 



JNg(z 



(e-**W[$(z)]u(l + 0{N- V ) ± eT^/[$(z)] 12 (l + ©(iV" 1 )) 



([$(z)]n ± e N ^z r {<S>(z)} 12 + ©(TV- 1 )) 
tN ^^p~F 1 {z) {{pF{z)) r ± e N ^ ( 



Kf) 



F{z) J pF{z) - pa 



0(N~ 



(4.34) 
(4.35) 



The leading term is analytic on £?; To check this, it is most convenient to use the expression 
in the first line using [e N ^/ 2 ] + = [c^^/ 2 ]- on B and the jump relation, ([$(^)]n, [$(z)]i2) + = 
{-z^{z)] l2 ,z-^{z)] xl )_. 
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Near j3 (and f3): Inside one can similarly obtain the strong asymptotics (1.16) of P, 1: m 
from j4°° (4.25) as in the previous calculations. We leave the (straightforward) calculation to the 
interested readers, mentioning only (1) that the relation = a/3/p (which is obtained from (2.12)) 
can be useful and (2) that we defined Cj at (1.15) such that 

y(z) = 2Cp(z - W' 2 (1 + 0(z - ~0)) . (4.36) 

In (1.16) of Theorem 1.1 we express the asymptotics in the scaled coordinate £(z) ~ (C^V) 2 / 3 (z— j3) 
taking only the leading linear term in the Taylor expansion of £(z) around (3. 



5 The post-critical situation: t > t c 



As t increases, (3 approaches to a as (3 ~ a + ac/t + 0(l/t 2 ) and b escapes to oo. We define A by 
the same equations (3.30) with the regions f2± as illustrated at the right side in Figure 8. The RHP 
for A is exactly the same as specified in eqs. 3.31 with the proviso that now T = B and hence the 
fourth condition in (3.31) is irrelevant. 

In this case, the solution (j> to the RHP at (4.1) is very simple. 



1 

-1 



z e Ext(B) , 
z G Int(B) . 



(5.1) 



Local parametrix at the double point: We define the local coordinate £(z) near f3 such that 

(z), z gB/j nlnt(B), 



C{z) = 2N- 



(z), zeD^nExtfB). 



(5.2) 



Under the mapping ( the contour B maps to the straight rays [0, e M / 4 oo) U [0, e~"/ 4 oo), and T + 
maps to the imaginary axis. See Figure 10 for the images of the contours under the conformal 
mapping £. 

where is a sufficiently small but fixed disc around z = (3 such that £(z) is one-to-one. We 
get the following expansion near z = f3. 



1 



CM _ 



(z-0)(l + O(z-0)), 7i 



W-o) 
a(b - (3) ■ 



(5.3) 



Inside we want to find V such that $(z) z ^ c ' 3 'P{z)z^' J3 satisfies the jump conditions of the 
RHP (3.31) for A; This leads to 



[V(z)} + = [P(z)l 



1 

e CW a /a 1 



z g r+ n : 



*0» 



V(z)=I + 0(l/N) 



-1 

1 



[V{z)l 



1 

-1 



z eBnBp, 

z G dDp . 



(5.4) 

(5.5) 
(5.6) 
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k 










Figure 10: (Left) Contour under ^-mapping; (Right) Plot of log |.F(C)| s -t- the brightness increases 
with the value. The dark dots are where the zeros of F are. (T is defined at (1.20).) 



A solution can be written by 
V(z) :- 

where 



H (z)F(((z)) , 

-1 

1 



F(C) 



1 







-1 

2iir 



H (z)F(((z)) 



ds 



1 

-1 



z e Ext(B) 0% 
z e Int(B) nD^, 



-too i> - C) 
1 



1 1 



/2VC 



1 




+ 0(C 



(5.7) 



(5.8) 



For any choice of H that is holomorphic and unimodular in D^, V defined by (5.7) satisfies the 
jump conditions in D^. We define the matrix Hq by 



H Q (z) := E(z) I 



' 


1 " 










) 



E(z) :=I 



1 liP r z~ 



V2ttN 



such that it satisfies, for z £ dU>p as N — > oo, 

Ho(z)F(C(z))E(z)- 1 =I + 0(azr 3 ) 



1 




(5.9) 



(5.10) 



The logic of how we choose E{z) is as follows: we require that i) S(z) is meromorphic in the extended 
complex plane with poles exactly at /? and oo; ii) the pole singularity of E(z) at f3 matches the pole 
singularity of 

I + 7SW,[<> <>}'■ (5 ' u) 

and iii) $(z)z~( r / 2 ) <T3 • E(z) • z^/ 2 ^ 3 has the same growth behavior as that of A{z) as z — > oo. 

This construction is based on the algorithm called "partial Schlesinger transform" introduced 
in [5] to compute the higher order corrections to the Riemann-Hilbert asymptotic analysis. In the 
critical case, we repeat the similar construction. 
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Now we define the strong asymptotics of A(z) by 



1 li(3 r z- r 


' 1 " 


) 


\/2ttN z-P 








A°°{z) := < 



-z-' r 




i J 3 


z r 


z 




z -i^-p(z)z^' T3 




1 7l /T 


' 


1 


V2tt7V z - P 








1 7l /T 





1 


V2tt7V z-P 









1 liP r z~ r 



V2ttN z ~P 






1 




z- r 


















-z r 



1+ 



<f>(z) z-i° 3 V(z) z l * 



z r 

z- r 

1 

-1 



z G Ext(B) \D /3 , 

z G Int(B) \D^, 
z G D«. 



(5.12) 



Remark. Note that, unlike the pre-critical case, the outer parametrix is modified from $ by a left 
factor. As we will see soon, if we used the unmodified outer parametrix $ in the definition of A°° 
we would have got the same strong asymptotics for P n ^N in Ext(B) with a larger error term, namely 
0(1/ vN). This error term turns out to be too big for the analysis in Int(23) because the leading 
behavior of P n ,N there is only of order 0(1/ 'y/N). 

Error analysis: We define the error matrix by 

£{z) :=A°°(z)A- 1 (z) . (5.13) 

One can show that the jump of £ is bounded such that the solution satisfies £ (z) = I + 0(A -3 / 2 ). 
Let us perform the computation on the contour, dDp H Ext(B). Below, + side is to the side of D^j. 

£{z) + (£{z)^ 

1 liP r z~ 



A 00 (z) + (A 00 (z)_)- 1 

$(z)z-2^p( z ) z §^ z -3 CT 3 



(using (5.7)) 
(using (5.10)) 



= z^H (z)F(az))[I- 



1 liP r z~ 



\/2ttN z-P 
I + (n- 3/2 ) , for z G dBp n Ext(B). 



z - 


-p 


"0 


r 










) 



"0 


1" 










) 


z~ 


5CT3 



z i a *<j>{z)- 1 (5.14) 

(5.15) 

(5.16) 
(5.17) 



A similar calculation gives the same error bound on Int(S) n <9B C . The jump matrices of £ on 
the other contours are exponentially close to the identity in large N and, therefore, £ satisfies a 
small-norm RHP and we conclude that 



A{z) = (l+0 (A^ 3/2 )) A°°{z) , 



(5.18) 



uniformly in the whole complex plane. 

As in the pre-critical case, the analysis of the asymptotic for the orthogonal polynomial can now 
be read off from our approximation in (5.12). 
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Away from B: The analysis of the asymptotic for the orthogonal polynomial can now be read 
off our approximation using that P n {z) = \Y{z)]u = Ane* 9 ' 2 ' outside 0±. Using the asymptotic 
behavior (5.18), we get the strong asymptotics outside the disk and Q± as follows (see also (3.27) 
and recall that tN = n — r). 



Pn(z) = 



I + O 



(AT 3 / 2 )) A°°(z) 



JNg(z) 



(l + 0(V 3 / 2 ))zV w ^), 
( liP r 



+ 



JNg(z 



z G Ext(S) , 
z G Int(B) , 



(5.19) 
(5.20) 



where I is given in (3.28). This holds uniformly over compact subsets of C \ B because T± can be 
moved as close to B as one wish. 



Near B but away from (3: In order to analyze the asymptotic on B we need to trace back the 
transformation that Y underwent within the region Q± . We thus obtain 



1.0- 






ff 05- 
i If ■ 








' \L c \ ' 




-i \ °: 




•"" \ 2 


^-0.5 - 












-1.0- 







Figure 11: The zeroes of P 59 with N = 30, a = 1, c = i, r = 0. Here t=f = 1.96 and t c w 1.8164. 
The thick line is where the zeroes converge. The solid (blue) line and the dotted lines are where 
(5.24) holds. As can be surmised from the figure, the dotted lines do not approaches B \ {f3} in the 
limit N —> oo; only the solid line approaches B \ {f3} and the roots seem to line up along that solid 
line. The shaded region indicates K. 



Pn(z) 



^A(z)e tN ^-i^ 



±- 



1 



z G fi±, 



(5.21) 



li 



34 



n = 55 




n = 70 




n = 60 




n = 80 



1.5 - 
/ f 0.5 - 




- 1 11 


o- 


Ml 2 


\\ 


-0.5 - 












-1.5 - 





77 
1.5 - 


- 65 




1.0 - 






1.5 1.0 10.5 
%j-0.5 - 

-1.0 - 

-1.5 - 


0.5 UT\ 


> 

1 .5 / 

/ 



n 

1.5- 


= 90 


/ 

/ f 0.5 - 

f i 

i 8 

i $ 




-ill 
\ \ -0.5 - 


Ml 2 


X l<0 - 




-1.5 - 




55,60,65, • • 


See the caption 



in Figure 11 for more details. 



where A = (I + 0(N 3 / 2 ))^4°°. A straightforward computation yields 



Pn{z) 



z r e tNg(z) , 1 _ 



z r e tNg(z) I e JV0( Z ) 



The leading term is analytic on B because [i g(.z)]± = [tg(z) + (f>(z)] T 



(5.22) 



Location of zeros: Here we explain the idea of the proof of Proposition 1.4 concentrating on 
the post-critical case. In order to locate the zeroes of P n we observe that they may only lie 
(asymptotically) where the following terms appearing above cancel each other: 



7i/? r 



1 



z G Q±. 



(5.23) 



Recall now that 7 = 0(s/N) and that ^t<p(z) < on both sides of B; this means that the zeroes of 
P n may possibly lie only within 57 + (i.e. slightly on the inside of B) on a region that lies within 
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= 0(lu.N/N). Therefore the location of the zeroes of P n (z) converges to the curve given 
implicitly by 



Re <t>(z) 



1 ,„ Ji? r 



InN 1 ~ a)(3 r 



(5.24) 



2N N y/2na(b- f3)z r (z- (3) 



Their (integrated) density is asymptotically determined by the difference of the imaginary parts of 
the two sides in (5.24) along the said curve. We remark that the above carries information on the 
location and density of zeroes only on parts of B that is away from the point /3. 

To prove the statement in Proposition 1.4 one starts by noticing that there is asymptotically no 
zero outside a strip of width 0(1/N) around the connected component of the set (5.24) that tends 
to B (at the rate 0{hxN /N)). This is seen by noticing that, where (5.24) is satisfied, the two leading 
terms in (5.22) for z G f2+ have the same order 0(N^ 1 / 2 ) and the error is of order 0(N~ 3 / 2 ). Next 
one uses the asymptotics of the polynomials together with some form of the argument principle, to 
show that the zeros actually lie in a strip of size 0(1/N) around the curve where (5.24) is satisfied. 
The Figures 11 and 12 clearly demonstrate such convergence. 

6 The critical case: t ~ t c 

For the critical case, tflg is the logarithmic potential for a signed measure and, therefore, 5J<^> may 
not satisfy the inequality 5i</> < near B. This does not cause any problem for us, because the 
region where the inequality is violated sits inside the domain of the local parametrix where we do 
not require the inequality. 

Following Definition 3.8 we define the (lens-opened) region f2± by the region enclosed by B and 
the steepest descent lines of di(j> from j3. We define A by the same formula as in (3.30) such that it 
satisfies the Riemann-Hilbert problem (3.31) except that there is no T \ B. 

Finitely away from b c = a + ^/c the jump matrices of A are all uniformly and exponentially close 
to the identity jump, save for the jump on B (3.31). We define $ by (5.1) to solve (4.1) as in the 
post-critical case. 

Local conformal coordinate: We remind the reader that <f> has been defined in Definition 3.4 
as a carefully defined anti-derivative of y(z) (which is defined in (2.30)), and possesses a jump 
discontinuity across B. As described in Section 2.3, the two critical points of 4> (fi fmd b) approach 
b c = a + y/c as N — > oo for the critical case. Below we intend to find a conformal map W which 
captures the important local behavior of 4> in a vicinity of b and (3. This is accomplished by defining 
a cubic polynomial 




(6.1) 



where 




(6.2) 



This cubic polynomial satisfies the following: 



Ps(Wp) = PsiWp) = Pai-Wp) = and 



8iP 3 {-W p ) = [<t>(b)] EM 



(6.3) 
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Lemma 6.1. For t near t c there exists a fixed disk D c centered at b c and the conformal map 
W : D c -> W(D C ) (? or eac fr *J satisfies 



8iP 3 (W(z)) = 



(z), zeD c nExt(6); 
0(z), zeD c nint(S). 



(6.4) 



The lemma is stated in Theorem 1 in [11] (with different notation: here we only note that the 
right hand side of (6.4) is an analytic function in D c because <p(z) has a jump discontinuity across 
B whereby <p + = —</>-). We also refer to Section 4.3 in [10] for another presentation. 

We define the scaled local coordinate £,(z), C,p and the scaled time coordinate s by 

1/6 7U2/3 

£(*):= JV 1/8 W(«), Cr-=N 1/3 Wp, s:^~AN 2 ^Wj^—— TW {t~t c ){l + 0{t-t c )). (6.5) 



such that we have 



2i^{zf + si{z) + \e^ 



'N<f>(z), zeD c nExt(6), 
-N<f>(z), z € D c nlnt(2?). 



(6.6) 



We will consider the vicinity t — t c ~ 0(N 2 / 3 ) such that s remains finite in the limit of large N. 
We define b* to be the root of W(z) = and a calculation (using 4>(b) + <j>([5) = 2</>(6*)) shows 

b* = b c + (t-t c )/(4b c )+(D((t-t c ) 2 ). 



Taking the derivative of (6.6) we have £'(z)\ z =b* = (—2is) 1 A r [y(6*)] Ext ( B ), which gives 



N l/3 



7 C 



(bc^C 1 / 6 



Later we need that, for z € O c and for a fixed N and s, we have 

^ = ^-bl){l + 0{z-bl))). 
See Figure 13 for how the various contours map under the mapping £. 



,1/3 



(6.7) 



(6. 



(6.9) 



Local RHP: We require that the local parametrix solve the RHP for A (3.31) exactly on D c ; i.e. 
we seek a matrix V(z) such that 

• ^(z)z~^ a3 V(z)z^ a3 solves the jump conditions of A within the disk O c ; 

• $(z)z"5' T 3-p( z ) z 5 CT 3 i s bounded in D c ; 
These conditions are satisfied if V satisfies 

-1 

1 



V(z) + 
V(z) + =V{z)_ 
V(z) + =V(z). 



V{z)- 



1 
-1 o j 

1 

e 2i(|c( 2 ) 3 + s «( 2 )+|Cl) i 

1 

e -2i(|C( Z ) 3 + S C(^)+|^) 1 



zeB 

zeT 



(6.10) 
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Figure 13: The contours in the local scaling ^-coordinate for t < t c (left) and for t > t c (right) 
under the mapping £. The branch cut (thick lines); the level lines of C 3 + S C + §Cj) from Qp 
and (a (dotted lines); the steepest ascent /descent lines (thin lines). 



The solution is not determined without specifying the boundary condition at z € dD c and, therefore, 
the solution can be written up to an undetermined holomorphic matrix H c (z) by 



V l l c ^ 3 iJ c (z)II(£(z))e I 3<> 3 , 



V(z) := 



1 

-1 



-1 

1 



z e Int(B) no c 
z e Ext(S) n D 



(6.11) 



Note that the jump conditions (6.10) for are satisfied for any holomorphic invertible matrix H c (z), 
if the matrix II solves the RHP with the following (in fact, only the first three) conditions. 



fn + (c) = n_(0, 
n+(c) = n_(c) 

n+(c) = n_(0 
[n(c) = / + 0(c ] 



1 o 

e 2i(|c 3 +sc) j 

i _ e - 2i (fc 3 +0 
o 1 



C -> oo; 



(6.12) 



See Figure 13 for the geometry of contours in ^-coordinate. The name of the contours in the local 
coordinate are given to match the corresponding contours under the mapping £; So ( G 8 implies 
that C £ £(£>) and so on. 



Remark. If we had not put the term H c (z) in (6.11) (i.e. H c (z) — I), we would have eventually 
found that 



V(z) = e -^ c ^ 3 U(^z)) e 1 ^ 3 =1 + 0(N- 1/3 ) 



z e 



(6.13) 
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This eventually determines the error of the strong asymptotics that turns out to be too large to see 
the leading effect in Int(B). This is why we introduce H c (z). 

The solution to the RHP (6.12) can be found in the literature about ordinary orthogonal polyno- 
mials [7] and is related to a particular solution called the Hastings-McLeod solution of the Painleve 
II equation. Especially one obtains 



n(C) = /- 



2iC 



u{s) q{s) 
-q(s) -u{s) 



O 



, u{s) = q'{sf - sq(s) 2 - q(s) 4 



(6.14) 



where q(s) is the so called Hastings-McLeod solution of the Painlev'e II equation: q"(s) = sq(s) 
2q(s) 3 . We note the asymptotic behavior: 



Q(s) = 



Ai(s) + O 



1 



' _i,3/ 2 
C 3 



3 l/4 



_ 2 3/2 
C 3 



1 + (V 3/2 n , 



<--s[l + 



(6.15) 



s — » —00; 



Strong asymptotics and error analysis: We follow the same construction of H c as that of H 
in the post-critical case. We define 



H c (z) = I 



1 



2i£(z) 



u(s) q(s) 
-q(s) -u(s) 



2N l l*{z-b* c ) \b 



u(s) q(s) 
-q(s) -u(s) 



(6.16) 



to exactly cancel the pole singularity of the first subleading term in (6.14). This is invertible for 
D c small enough. The conjugation by (z/b* c )i a3 is to later match the growth behavior of the outer 
parametric A°° at z = 00 that we define soon. Then we have 



where we have defined the following shorthand notation for convenience. 



z e 



S°°{z) 



u(s) 



2N 1 / 3 (z-b* c ) -q(s)(b*) r e 1 ^^ 
Now we define our strong asymptotics for A(z) by 
'<&(*) S°° (2), 
A°°(z) := { $(z)z"^ 3 



q(s)(b*)- r e- 1 -^ 
-u(s) 



(6.17) 



(6.18) 



1 
-1 

,$(2) z~^" 3 V{z) z%° 3 , 



I 



2N 1 / 3 (z-b*) 
t $(z) z-i a3 V(z) Z5" 3 



Z 5' T3 5 00 (Z)Z-5 <T3 

-u(s) 



-1 

1 



z^ 



-q(s)(b* 



q(s)(b*) r e 1 ^ 
e ''^ u(s) 



z e Int(B) \D C , 
zeExt(6)\D c , (6.19) 

z e D c , 

j*(*)> z$I 
z el 



(6.20) 
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As in the pre-critical case, we define the error matrix by 

£{z) :=A 00 (z)A- 1 (z). (6.21) 

Let us evaluate the jump of the error matrix on Int(£>) n <9B C . Below, + side is to the side of D c . 

£{z) + {£{z)_)^ = A™{z) + (A°°{z)_)^ 

= $(z)z-5 CT3 P(z)zi' T3 (5 00 (z))- 1 $(z)- 1 



(using (6.11)) 
(using (6.17)) 



1 
-1 

1 
-1 



i + o(n- 2,z ^ 



z-^e-n^H^n^z^S^z^iS^iz))- 1 

y-i^ z l^S' yo {z)z-^ a (i + O (a^ 2/3 )) z^ 3 (S°°{z)) 
for z e Int(S) n <9B C . 



-1 

1 



-1 

1 



(6.22) 



A similar calculation gives the same error bound on Ext(£>) (1 dD c . The jump matrices of £ on 
the other contours are exponentially close to the identity in large N and, therefore, £ satisfies a 
small-norm RHP and we conclude that 



A{z) =(l + (at- 2 / 3 )) A°°(z) , 



(6.23) 



uniformly in the whole complex plane. 

To find the strong asymptotics of the orthogonal polynomial we need to undo the transformations 
that lead us from Y(z) to A(z) and then rewrite it in terms of the approximation A°°(z). We skip 
the calculations as they are similar to the ones for the pre-critical and the post-critical cases. We 
note that the final form of the asymptotics contain b* and 7* which can subsequently be replaced 
by b c and j c using the approximations (6.7) and (6.8). This gives Theorem 1.3. 

To show the Proposition 1.5 near j3 we also need to evaluate P n ,N inside D c , that involves the 
Painleve II parametric II. We only need to use that 11(C) m (6-12) is bounded uniformly over £ in 
any compact subsets of C. 



7 Proof of Theorem 1.6 

First we relate the norming constant of P n .N with respect to the weight w n ^(z) dz with the norming 
constant, h n , with respect to e -Ar< ^' dA(z). Let us define the former by (recall the definition of Vj 
in (3.13)) 

h n := <p P n N {z)w n , N Az = —— ; ~ — — . (7.1) 

Jr ' det[i^ +:) 'Jo<ij<n-i 

Proposition 7.1. We have (note n = Nt + r) 
T(Nc + n+l) h n 

2lN Nc + n+l Pn+1N (Q) V- Z > 

\ \NJJ V ' \ 2N Pu+un(0) V ' 
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The proof can be found in Appendix E. 

The importance of Proposition 7.1 relies in the fact that we can find the asymptotics for h n in 
terms of the RHP already set up. Indeed it follows from eq. (3.15) that 



h n = -2m lim z n+L [Y(z)] 12 (7.4) 

z— f oo 

For the post-critical case, using (5.12), we have 

1 /2tt 



h n = -2m \imz n+1 [Y(z)} 12 = (/ + O (n^)) i^ 7l/ Te^, (7.5) 
and using g(0) — £ (3.27) (because G Int(B)) we have 

*W(0) = (l + (N-^)) ^/Te^") = (l+0 (N-^)) ^L 7l /Te^. (7.6) 
This gives h n / P n+ \^{Q) rs — 27ri and, from Proposition 7.1 we obtain 

h n = (l+0(7V- 3 / 2 ) ) 2*J^j(c + t )r+l/2 e Ar[( t + C )ln(t +C )-(t+ C )]_ (7 ?) 

Understanding that p = \/c + t and ^2D = (i + c ) ln(t + c) — (t + c) for the post-critical (and critical) 
case, we obtain Theorem 1.6 for the post-critical case. The calculation for the critical case goes 
exactly the same except that the error bound gets worse to 0(N~ 2 / 3 ). 

Let us evaluate h n for the pre-critical case. Using A°° in (4.25), we obtain 



h n = —2iir lim z 



n+l 

z— >oo 



e tN^a 3 



-2m 



(i + O (7V- 3/2 /(l + |z|))) A°°(z) e tN ^)-D^ ^ (7.8) 

i{ap) r ^- p e NU (l + o[N- z ' 2 )y (7.9) 

The z-dependence in the error O (iV -3 / 2 / (1 + \z\j) comes because the non-trivial jumps of the error 
matrix £ (5.13) are on a bounded set of contours. 

Using the identities (that follow from (2.17) and the last equation in (A. 2)) 



P F(0) = \0\ = ap + k/ P =(c + t)a/p, 

p f'(o) = (ifl - mum) = */{a 2 p + « 



pn 1 (7.10) 

OL z C + t 



we have 



P„ + i(0) = (1 +0(l/N))-fiF 7 (fi)(i,F(»)r"e m ' m 

Mi +0 (iw)^E (c+t) ;;r' e^°'. (7 - u) 

V or c + t p r+L 

Collecting the above results with Proposition 7.1 we obtain the following asymptotic result. 

*- ( 1+ K^)) i(c+tr ^^ e " [< ' +c>M ' + ^ M CT (7 ' 12) 

= (l + O 2n \p^~ P 2r+1 e N{tl ~ t9 ^ ) ~ { - t+c)+{t+c)Xo&( - t+c)) . (7.13) 
We further need the following result. 
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Lemma 7.2. Defining £20 ■— 2tiftg(z) — Q(z) for z £ dK. Then we have 

izo = t(t - n 9 (0)) -(c + t) + {c + t) log(c + t). (7.14) 

Proof. We take the limit z — > at the equation in Lemma 2.12. As z — > we have Shaded) = 
(c + t)/z + 0(l) -> 00 and S hack (z) z = c + 1 + O (z) . From W OP (2) = 3lV(z) - 2tUg(z) +U we have 

2U OP {z) = -2cloga + 2(c + t)log|z| - 4t%(0) +2« + C(z), (7.15) 
+W(5 back (2)) - |S back (z) - z| 2 = -2c log a - 2i»ff(0) + 2^ 2D (7.16) 
- 2(c + t) log((c + t)/\z\) + 2(c + 1) + 0(2), (7.17) 

as z — > 0. Equating the above two, we obtain the proposed equation in the limit z — > 0. □ 

Combining (7.13) with Lemma 7.2 we obtain Theorem 1.6 for the pre-critical case. 



A Proof of Lemma 2.5 

Using the deck transformation (2.11), the injectivity of / on |u| > 1 is equivalent to the condition: 
deck(D ) n D =0. Since the deck transformation is a Mobius transform and maps the unit circle 
to the circle that passes through deck(±l), the above injectivity is equivalent to the condition: 
— 1 < dcck(±l) < 1, which gives the equation (2.10). 

Suppose that / : D — > K c is a conformal mapping. (In the last step of the proof, we will check 
(2.10 ) to confirm this.) Since f(v) G dK when v € <9D it follows that S{z) = S{f(v)) = f{l/v) = 
f(v) = z over z £ dK. So the item i) in Definition 2.1 is satisfied. To satisfiy ii) and iii) we note 
that 

f(l/v) —¥ oo only when w->l/aorw-4 0. (A.l) 

It means that the full (multi-valued) analytic continuation of S(z) over the Riemann surface has 
two poles at /(1/a) and /(0) = 0. To satisfy ii) in Definition 2.1, since a ^ 0, we have f(l/a) = a. 
Therefore, the conditions ii) and iii) in Definition 2.1 are equivalent to the following three equations 
for the the three unknowns p, k, a: 

, , . pa 2 - p + K 
/(V") = 



a(a 2 — 1) 

res S(z) dz = res f(l/v) f(v)dv =^( P + n = c, ( A - 2 ) 

z=a v=l/a OL z \ \\ — U A y) 

- res S(z) dz = - res f(l/v) f(v)dv =^+p 2 = c + t. 

z— oo v—oo (X 

We show that these equations have a unique solution within the proposed parameter space. By 
eliminating p and k from equations (A.2) we obtain 

t + a 2 a 2 (l~a 2 )(t~a 2 a 2 ) 

P=—~ > K =~ £ A -3 

2a a 2a a 

and a must satisfy 

* <A,, 
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This is a third order equation in X and its discriminant in X is given by i 2 A/(4a 10 ) where 

A := (a 2 + 4c + 2i) 3 - 27tV. (A.5) 

This is strictly positive for any c > because d c A = 6(a 2 + 4c + 2i) 2 > and A| c= o = (a 2 — 
i) 2 (a 2 + 8i) > 0. Therefore the equation (A. 4) has three real roots among which exactly two of 
them arc positive, because the sum of the three roots is positive and the product of the three roots 
is negative. For < t < t c = a(a + 2y/c) we have c > (a 2 — t) 2 / (4a 2 ), and 

P(0) = t 2 /{2a 4 ) > 0; P(l) = -2 [c - (a 2 - t) 2 /{Aa 2 )} /a 2 < 0. (A.6) 

So, there is exactly one real root of (A. 4) at < X < 1 for < t < t c . We take a = yX € (0, 1) 
as the unique solution of (A. 2). Also by P(t/a 2 ) — —2ct 2 /a 6 < we have a 2 < t/a 2 . Therefore k 
in (A. 3) is positive for our choice of a. 

In summary the system of equations (A. 2) has exactly one solution (p, k, a) in K + x R + x (0, 1). 
Lastly, we need to check whether our solution satisfies (2.10). From (A. 3) it follows that 

k , o. t — a 2 a 2 n , . 

- = (l-a 2 )— — <l~a 2 . (A.7) 
p t + a z a z 

To conclude, the region K defined through the conformal map / has the Schwarz function that 
satisfies Definition 2.1. 

B Proof of Lemma 2.9 

Proof. Using the definition (2.21), S (z) + S\, a( ±(z) defines a rational function because it is continuous 
on B; Crossing B the boundary values satisfy [F(z)]± = [dcckoF(z)] T by Lemma 2.8, and therefore 
[>S'( 2: )]± = [Shack (^)] =p- To determine the rational function we only need to identify pole singularities. 

- S( z ) = f(l/F(zj) has a pole at p only if i) F(p) = 1/a or ii) F(p) = 0. 

- Sbadi(z) = f (1 / ^decko F (z)) has a pole at p only if iii) deckoF(p) = l/aoriv) deckoF^) = 0. 

- i) and iii) are complementary, and p = /(1/a) for any of them. 

- ii) and iv) are complementary, p — /(0) = for any of them. 

- Since we know that S(z) has a pole at a with residue c, i) must hold with p = a. 

- ii) is impossible because, F(0) — \/3\/p ^ by using the fact that is to the right of the 
branch cut B. 

- So iv) holds, i.e. deck o F(0) — 0. The residue of the pole of Shack at z = is given by 

res Sback(z) dz = res f(l/v) df(v) = — res f(l/u)df(u) = — res S(z) dz = t + c. (B.l) 

z— v— u— 00 z=oo 

- Lastly, S b ack(oo) = /(1/deck o F(oo)) = /(l/deck(oo)) = /(1/a) = a. 
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Above items uniquely determine 

S(z) + S hlM:k (z) = a + c— + — . (B.2) 
z — a z 

The difference satisties 

S(z) — S\, ack (z) = (F( z ) — decko F(z)) ■ (rational map symmetric under F(z) -H> decko F(z)) 



(z — 0){z — {$) ■ (rational function). (B.3) 

Using the pole structure of S and <Sb ac k at a, 0, oo, the rational function is determined such that 
S{z) — S'back^) = ~y( z ) has the pole structure as proposed. □ 

C Proof of Lemma 2.12 and Lemma 3.3 

The proof is based on the anti-holomorphic involution symmetry of v 1 /v. Below, we use the 
definition S hlick (f (v)) = /(l/deck(u)) (2.21), and set v = F(z). 

y(w)dw= f (Sback(O-S(O)dC 



deck(t>) — > w; 
Note deck(t> / 3) = vp 



/(l/deck(u)) d/(deck(«)) + / /(1/u) df(u) 

d.eck(v) pvp 

f(l/w)df(w) + f(l/u)df(u) 

deck(if) 

f(l/u)df(u). (C.l) 



Choosing any point vq € dD the last equation can be written by integration by parts: 

y(w)dw= / + / /(l/«)d/(«) (C.2) 



/3 



v rl/vo 

/(l/«)d/(u)+ / /(l/u)d/(«) 

Jl/deck(u) 



dcck(v) 

f(l/v)f(v) . (C.3) 



Using l/« = wo and /(l/deck(i>)) = 5 ba ck(^), we have (defining z := f(v )) 

S(z)dz + zS h nck(z) + zS hack (z) -2|z | 5 



2W OP (z) = 29? ^ y(s) ds = 25R ^ 



-20 /"^O 



Shack 0) 



= «(*) + W(S back (z)) - (Z ~ S h , ck (z))(z - S back (z)). (C.4) 

This proves Lemma 2.12. 
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Below, we present the proof of Lemma 3.3. We remark that, if v = F(z) satisfies deck(w) = 1/v 
then the left hand side of (C.2) is purely real, which can be confirmed by 



mMf(v) 



(C.5) 



pVo rl/vo r -| deck(u) 

y(w)dw = /(l/u)d/(u)+ / f(l/u)df(u) 

13 Jv Jl/dcckO) 

/(l/u)d/(u)+ / f(l/u)df(u) + f(v)f(deck(v))-f(vS)f(v ) (C.6) 



-23? / S(z) dz + \z\ 2 - \z \ 2 = U(z). 

J z 



(C.7) 



Since b satisfies deck(F(6)) = 1/F(b), we have <j>(b) =U(b). Since a point on the steepest lines from 
P has purely real <f> value, b must be on a steepest line from /3. 

D Proof of Proposition 3.2 



Up to a sign due to a permutation of columns we have 
det[^ +J -] < iJ < n _ 1 = (-lr^-^det 



i_j (z - a) Nc c- Naz dz 

z Nc+l 



(D.l) 



0<i,j<n— 1 



All the determinants in this section is for n by n matrix, i.e. < i,j < n — 1. Due to Lemma 3.1 
the latter determinant is given by 



det 



l _ j (z - a) Nc e- Naz dz 
z N ^ 



= det 



JJ z\z - af\z - a\™ D er N W* &A{z) 



n-l 2tN Nc+k+l 

W r(Nc + k + 1)' 

(D.2) 



J fc=0 



Finally, the determinant on the right-hand side is strictly positive because 



det 



z-a) 3 \z-a\ 2Nc e~ N ^' 



dA(z) 



dot 



if 



z'z-* |z — a 



2JV C „-JV|*|' 



dA(js) 



> (D.3) 



where the equality is due to the fact that the columns of the two matrices are related by a unimodular 
triangular matrix, whereas the inequality follows from the positivity of the measure. Finally, since 
r(z) has no zeroes (and no poles since Nc + k + 1 > 1), the nonvanishing follows from (D.2). 



E Proof of Proposition 7.1 

From the determinantal expressions for P n ,Ni defining /i^ := $ z % z?\z — a| 2Arc e _Ar ' z ' d 2 z, we have 

det[fj.ij]o<i,j< n 



det \pij]o<ij< n -i 

The denominator already appeared in the proof of Proposition 3.2 and is given by 
det[^]o<ij<n-i = (-l) n(n - 1)/2 det[i/ i+j ] <i,j<n-i \[ ^^^^ 



fc=0 



(E.l) 



(E.2) 
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For the numerator in h n we have 

(, \Nc p -Naz 



det[^ +J ]o<ij< n = det 



yNc+l 



dz 



( _ 1) n(„+l)/2 TT ^ 



z 4+1 (^^) fc |z - a| aJVc e- w l«l a dA(*) 



0<i,k<n 



= ( _ 1) n(«+l)/2 TT 

£=0 v ; 



= (-1 



D" ( " +1)/2 n fWTTTT) ("ir 1 W(o) det 



[Mifc]o<i,fe<n 



(E.3) 



Dividing both sides by det[j/; + j]o<ij< n -i and using eq. (E.2) we prove the first equality in Propo- 
sition 7.1. 

deth +J ]o<^< w 2iJV^+"+ 1 
det[u i+ j\o<i,j<n-i l(Nc + n + l) 



K = 



(E.4) 



The last equality comes from the Stirling approximation formula 



T(z) = ] / T z*e-*(l + 0(z-i)). 

We find (recall that n = tN + r) 

r(N(c + t)+r + l) f (1\\ ( . +1 T(N(c + t)) 

NN(c+t)+r+l I \N ] J ' ' N N i c + l ) 



(E.5) 



( 1 + (±y\ ( C + ty ] / 2 ^±^e N M*+c)-(*+o] (E 6) 
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